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Abstract 

Holographic  Interferometry  and  computer- 
assisted  tomography  (CAT)  are  used  to  determine 
the  transonic  velocity  field  of  a  model  rotor 
blade  In  hover.  A  pulsed  ruby  laser  recorded 
<10  lnterferograms  with  a  2-ft-dlam  view  field  near 
the  model  rotor-blade  tip  operating  at  a  tip  Mach 
number  of  0.90.  After  digitizing  the  lnterfero¬ 
grams  and  extracting  fringe-order  functions,  the 
data  are  transferred  to  a  CAT  code.-  The  CAT  code 
then  calculates  the  perturbation  velocity  In  sev¬ 
eral  planes  above  the  blade  surface.  The  values 
from  the  holography-CAT  method  compare  favorably 
with  previously  obtained  numerical  computations  In 
most  locations  near  the  blade  tip.  The  results 
demonstrate  the  technique's  potential  for  three- 
dimensional  transonic  rotor  flow  studies. 


Notation 

a0  -  3peed  of  sound,  ft/sec 

A  -  wave  amplitude 

C  *  blade  chord,  ft 

f(x')  -  filter  function 

I  -  lrradlance 

k  -  Gladstone-Dale  constant,  ft^/slug 

L  *  path  length,  ft 

n  •  refractive  Index 

n0  -  ambient  refractive  Index 
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N  *  fringe-order  number 

R  *  spanwlse  coordinate,  ft 

R0  -  blade  span,  ft 

t  -  hologram  amplitude  transmittance 

Uc  *  reconstruction  wave  complex  amplitude 

Uj  *  transmitted  wave  complex  amplitude 
Up  -  reference  wave  complex  amplitude 
UQ1  -  ambient  object  wave  complex  amplitude 
U02  *  test  object  wave  complex  amplitude 

Y  -  perturbation  velocity,  ft/sec 

x'  -  projection  coordinate,  ft 

X  -  chordwlse  coordinate,  ft 

Y  -  height  above  blade  centerline,  ft 

B  -  film  proportionality  constant 

Y  -  ratio  of  specific  heals 

*  optical  path-length  difference  (OPD) 

9  *  field  projection  angle,  deg 

X  -  laser  wavelength,  ft 

2 

p  *  air  density,  slug/ft 

pQ  *  ambient  air  density,  slugs/rt^ 

(2  *  blade  rotational  speed,  rpm 


Introduction 

On  many  helicopters,  the  rotor  blade's 
advancing  tip  encounters  transonic  flow  during 
forward  flight.  St  these  high  Mach  numbers,  the 
rotor  blade's  performance  suffers  from  compressi¬ 
bility  effects  that  often  cause  shock  waves  to 
form  near  the  blade  tip;  the  shocks  can  extend  to 
the  acoustic  f ai — field.  Through  theoretical  and 
computational  investigations,  researchers  attempt 
to  understand  the  local  shock  generation  of  high- 
tip-speed  rotors  and  its  propagation  to  the  far- 
field.  However,  because  of  the  problem's  complex¬ 
ity  and  the  difficulty  of  obtaining  detailed 
experimental  information  about  the  flow,  accurate 
means  for  confirming  transonic  rotor-blade  designs 
have  been  notably  lacking. 

Shock  waves  cause  a  number  of  aerodynamic, 
dynamic,  and  acoustic  problems  on  high-speed  hells 
copter  rotor  blades.  First,  the  shock  rapidly 
increases  the  aerodynamic  drag  through  energy 
dissipation,  flow  separation,  and  wave  effects. 
Second,  local  shocks  cause  sudden  large  changes  in 
pitching  moment  which  can  excite  various  blade 
torsional  modes.  As  the  blade  rotates  in  forward 
flight,  its  Mach  number  and  angle  of _ attack  vary. 
The  shock  appears  on  the  advancing  side  of  the 
rotor  disk  and  often  results  in  large  chordwise 
movements;  these  movements  can  be  in  opposite 
directions  on  the  upper  and  lower  surfaces  as  the 
Mach  number  and  angle  of  attack  change  during  each 
revolution.  The  changing  shock  positions  on  the 
upper  and  lower  surfaces  cause  an  unsteady  load¬ 
ing,  which  produces  fluctuating  pitching  moments. 
These  moments  can  cause  unexpected  blade  motions, 
oscillating  loads  on  pitch  links,  and  vibrations 
throughout  the  entire  aircraft.  Third,  shock 
waves  on  an  advanclng-blade  surface  can  "delocal¬ 
ize"'  and  extend  directly  to  the  far-field.  Large 
amounts  of  acoustic  energy  radiate  in  front  of  the 
helicopter  near  the  tip-path  plane.  This  heli¬ 
copter  impulsive  noise  is  annoying  in  general  and 
too  easily  detected  in  military  applications. 

In  an  attempt  to  describe  the  transonic  rotor 
flow  field  and  to  resolve  the  problems  associated 
with  it,  promising  theoretical  models2  and  numeri¬ 
cal  codes^’11  have  been  developed.  The  numerical 
codes  compare  favorably  with  blade-pressure  mea¬ 
surements,"’  but  are  not  yet  verified  at  points 
away  from  the  blade's  surface. 

Previous  attempts  to  measure  the  flow  field 
have  been  severely  limited.  Pressure-instrumented 
alrfoll3  are  expensive  and  difficult  to  fabri- 
cate- -espec i a  1 1 y  in  scale  models.  In  addition, 
pressures  can  be  measured  only  at  the  blade  sur¬ 
face.  Hot.-wite  anemometry  requires  that  a  probe 
be  positioned  within  the  field,  therefore  disturb¬ 
ing  the  flow.  Laser  velocimetry  requires 


flow-seeding  and,  when  shock  waves  are  present,  it 
is  uncertain  whether  the  seeds  follow  the  flow 
faithfully.  Both  hot-wire  anemometry  and  laser 
velocimetry  can  take  only  point-by-point  measure¬ 
ments,  requiring  large  amounts  of  running  time  to 
survey  the  rotor's  three-dimensional  field,  a 
distinct  disadvantage  for  rotor  testing. 

Schlieren  and  shadowgraph  photography  provide  only 
a  qualitative  two-dimensional  representation  of  a 
three-dimensional  flow.  And  Mach-Zehnder  inter¬ 
ferometry  provides  quantitative  information,  but 
is  extremely  difficult  to  use  in  a  large-scale 
experiment.  Clearly,  another  experimental  tech¬ 
nique  which  overcomes  these  limitations  must  be 
employed . 

Holographic  interferometry  is  an  effective 
diagnostic  technique  for  making  transonic  flow 
measurements.®  Previous  investigations^’®  in 
which  two-dimensional  flow3  over  airfoils  were 
studied  show  that  accurate  quantitative  informa¬ 
tion  is  obtainable  using  holographic  interferom¬ 
etry.  However,  the  transonic  flow  around  a  heli¬ 
copter  rotor  blade  is  three-dimensional  and 
requires  a  tomographic  technique  to  compute  the 
correct  flow  information  from  several  interfero- 
grams.  To  date,  most  applications  of  this  tech¬ 
nique  have  been  limited  to  axisymmetric  flow  or  to 
simple  three-dimensional  flow  with  a  small  model 
under  ideal  laboratory  conditions.0,10 

This  paper  discusses  the  procedures  necessary 
to  obtain  quantitative  measurements  of  a  tran¬ 
sonic,  three-dimensional  flow  field  near  a  rotor- 
blade  tip,  using  holographic  interferometric  data 
and  CAT.  Though  most  helicopter  rotor  problems 
caused  by  shock  waves  occur  during  forward  flight, 
this  experiment  investigates  the  steady  problem 
(hover),  which  simulates  many  physical  phenomena 
of  forward  flight.'1  The  method  for  recording 
Interferograms  and  example  lnterferograms  is 
Included,  and  the  steps  that  must  be  followed  in 
extracting  quantitative  information  from  the 
lnterferograms  are  outlined.  The  technique's 
potential  for  measuring  three-dimensional  tran¬ 
sonic  rotor  riows  is  demonstrated,  and  the  results 
it  yields  are  compared  with  those  from  previously 
performed  numerical  computations. 

Background  Concepts 

For  the  experiment  to  be  successful,  it  is 
necessary  to  1)  record  high-quality  i nt erferograms 
near  a  rotor-blade  tip  from  multiple  viewing 
angles,  and  2)  implement  a  suitable  CAT  code  with 
the  i nterferogram  data.  Familiarity  with  holog¬ 
raphy,  holographic  interferometry,  and  computer- 
assisted  tomography  principles  provides  the  neces¬ 
sary  1  ns i gilt  for  understanding  this  technique. 
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Holography 

Holography  13  a  two-step  Imaging  process  In 
which  diffracted  light  waves  are  recorded  and 
reconstructed.  ’  J  The  first  step  is  recording, 
or  storing,  the  hologram.  This  is  accomplished  by 
dividing  a  single  coherent  laser  beam  Into  two 
beams  and  exposing  a  photographic  film  to  the  two 
light  waves,  as  shown  in  Fig.  la.  The  object 
wave,  which  Is  the  wave  containing  the  flow  infor¬ 
mation,  passes  through  the.  measured  field  (the  air 
near  the  blade  tip  In  this  experiment).  The  sec¬ 
ond  wave,  the  reference  wave,  passes  around  the 
field.  By  adding  the  coherent  reference  wave  to 
the  object  wave,  the  photographic  film  records  a 
high-frequency  interference  pattern.  Once  the 
film  is  developed,  it  is  known  as  a  hologram, 
which  i 3  a  complicated  diffraction  grating. 

The  second  step  in  holography  is  reconstruc-  . 
tion,  or  playing  back  the  hologram.  This  is 
accomplished  in  two  ways.  First,  a  reconstruction 
wave  identical  to  the  reference  wave  illuminates 
the  hologram  (Fig.  1b).  The  hologram  diffracts 
the  reconstruction  wave  and  produces  a  replica  of 
the  original  object  wave,  forming  the  original 
object's  virtual  image.  In  the  second  method  of 
reconstruct  ton ,  a  reconstruction  wave  conjugate  to 
the  reference  wave  Illuminates  the  hologram 
(Fig.  1c).  The  hologram  dlfTracts  the  conjugate 
wave  forming  the  original  object's  real  image. 

The  real  image  may  be  photographed  without  the  use 
of  a  lens  by  placing  a  sheet  or  photographic  film 
in  the  real  image  space. 

Several  Important  characteristics  of  hologra¬ 
phy  are  applicable  to  the  experiment  at  hand. 

There  are  very  few  geometrical  constraints  In  a 
^olographic  optical  system;  thus,  holography  can 
be  applied  In  a  large-scale,  nonlaboratory  envi¬ 
ronment.  Note  that  recording  and  reconstruction 
of  the  hologram  can  be  done  in  different  locations 
if  the  reference  wave  i3  reproducible.  Thl3 
allows  the  reconstruction  to  be  done  In  a  labora¬ 
tory,  far  from  the  harsh  environment  in  which  the 
hologram  was  previously  recorded.  The  reference 
wave  serves  only  as  a  method  of  recording  and 
reconstructing  the  object  wave.  Thus,  a  hologram 
does  not  produce  quantitative  information  about 
the  field  of  interest.  To  obtain  quantitative 
information  (in  the  form  of  interference  fringes) 
an  lnterferogram  must  be  recorded. 


ferogram  is  recorded  by  first  exposing  a  photo¬ 
graphic  film  to  i.  reference  wave  and  an  "undis¬ 
turbed"  object  wave.  Later  in  time,  the  same 
photographic  film  is  exposed  to  a  reference  wave 
and  to  a  second  "disturbed”  object  wave. 

When  the  holographic  lnterferogram  is  recon¬ 
structed,  the  virtual  or  real  image  shows  the 
object  (the  transparent  field)  with  an  interfer¬ 
ence  fringe  pattern.  The  fringe  pattern  repre¬ 
sents  the  difference  between  the  "undisturbed"  and 
"disturbed"  flow  states.  The  lrradiance  of  the 
reconstructed  wave  Is  proportional  to 


which  can  be  written  as12’1^ 

I  -  2A2[1  *  cos(A*)]  (1) 

Equation  (1)  represents  the  lnterferogram  with  a 
fringe  pattern  of  dark  and  bright  fringes  of  con¬ 
stant  optical  path-length  difference  (OPD)  A<|>, 
where  A$  is  given  by 

A*  -  J [n(x, y,z)  -  nQ]  -  NA  (2) 

To  determine  the  flow-field  properties,  the  line 
Integral  of  Eq.  (2)  must  be  Inverted  and  solved 
for  n(x,y,z),  the  refractive  index  at  a  specific 
point  in  the  field. 

In  a  two-dimensional  flow  (l.e.,  the  flow 
over  a  fixed  airfoil  in  a  wind  tunnel),  the  evalu¬ 
ation  of  Eq.  (2)  is  simplified,  since  the  refrac¬ 
tive  index  is  constant  across  the  width  of  the 
test  section  L,  Eq.  (2)  becomes 

n(x,y)  "  n0  *  NA/L  (3) 

In  a  two-dimensional  flow,  the  fringes  on  an 
lnterferogram  are  contours  of  constant  refractive 
index  and  the  refractive  index  at  any  point  in  the 
field  can  be  determined  from  a  single  interfero- 
gram.  However,  3inre  the  transonic  flow  near  a 
rotor-blade  tip  is  three-dimensional,  Eq .  (3) 
cannot  be  used.  To  Invert  Eq.  (2)  and  solve  for 
the  refractive  index  at  a  specific  point  in  the 
field,  computer-assisted  tomography  (CAT)  must  be 
employed . 

Computer-Assisted  Tomography 


Holographic  Interferometry 

Holographic  Interferometry  is  the  interfero¬ 
metric  comparison  of  two  object  waves  recorded 
holographies! 1 y  (see  the  Appendix  for  further 
detail).  In  this  experiment,  the  two  object  waves 
are  recorded  sequentially  in  time  with  double¬ 
exposure  holographic  Interferometry.  The  Inter- 


Tomography  is  a  mathematical  technique  for 
reconstructing  a  three-dimensional  field  from  ita 
two-dimensional  projections  (see  References  1 H 
and  15  for  a  wide  variety  of  applications).  A 
projection  of  a  three-dimensional  field  is  the 
fringe  pattern  recorded  on  an  lnterferogram.  All 
methods  require  multidirectional  projection  data 
of  the  field.  Tomographic  codes  develop  in  two 


directions:  1)  iterative  algebraic  reconstruction 

techniques,1®  and  2)  Fourier  transform  techniques. 
A  version  of  the  latter  method,  termed  filtered 
back-projection,1^’1®  appears  most  suitable  for 
this  application. 

Most  Fourier  transform  techniques  employ 
back-projection.  Projection  data  from  the  field 
are  recorded  in  one  plane  at  several  azimuthal 
angles  around  the  field.  For  example,  one  projec¬ 
tion  of  a  uniform  absorbing  disk  is  shown  in 
Fig.  2a  (taken  from  Reference  19).  Beyond  the 
disk  boundary  (no  path  length  through  the  disk), 
the  light  ray's  OPD  is  unchanged,  producing  no 
interference  fringes.  Near  the  disk  boundary 
(short  path  length  through  the  disk),  the  light 
ray's  OPD  is  changed  slightly,  producing  a  few 
interference  fringes.  And  near  the  disk  center 
(long  path  length  through  the  absorbing  disk),  the 
light  ray's  OPD  is  changed  substantially,  produo- ‘ 
i  ng  several  interference  fringes.  'limilar  projec¬ 
tions  (fringe  number  vs  position)  at  different 
azimuthal  angles  are  also  recorded.  Faoh  projec¬ 
tion  Is  then  back-projected,  or  smeared  back  along 
the  direction  In  which  it  was  recorded  (Fig.  2b). 
Values  are  added  point  by  point  to  form  a  recon¬ 
struction  of  the  field.  Unfortunately ,  simple 
hack-projection  produces  an  undesirable  spoke 
pattern  which  severely  degrades  the  quality  of  the 
reconstructed  field. 

To  eliminate  the  spoke  pattern,  the  pro¬ 
jections  are  filtered.  A  one-dimensional  con¬ 
volution  (indicated  by  an  asterisk)  is  performed 
between  each  projection  and  an  appropriate  filter 
function  (see  Reference  20  for  a  discussion  of 
filter  functions)  before  back-projection,  as  shown 
in  Fig.  ja  (taken  from  Reference  19).  Each  fil¬ 
tered  projection  is  then  back-projected  over  the 
reconstruction  space  (see  Fig.  3b).  The  negative 
side-lobes  introduced  by  the  filter  eliminate  the 
spoke  pattern  during  the  point-by-point  addition 
process.  With  many  projections,  this  technique 
yhields  an  accurate  reconstruction  of  the  original 
field. 


Procedure 

Several  steps  must  be  performed  to  quantita¬ 
tively  reconstruct  the  three-dimensional  transonic 
field  near  a  model  helicopter  blade  tip.  First, 
several  holographic  i nt.erf erograms  must  be 
recorded  along  planes  perpendicular  to  the  rotor 
tip-path  plane  at  various  azimuthal  angles 
(Fig.  A).  Data  must  then  be  extracted  from  the 
i nterferograms .  This  can  be  done  1)  manually, 

2)  by  U3ing  a  graphic  tracing  tahlet.,21  or  3)  j,y 
using  a  system  that  digitizes  the  1 nterferograms 
and  extracts  fringe-order  numbers.  The  digital 
I nterferogram  evaluation  technique  was  used;  it 


Op 

will  be  presented  in  detail  by  Becker  and  Yu. 
Finally,  the  data  are  transferred  as  input  to  a 
tomography  code,  which  computes  the  refractive 
index  at  specific  points  in  a  horiz.ontal  plane 
above  the  blade  surface.  This  procedure  is 
repeated  in  several  planes  to  yield  a  reconstruc¬ 
tion  of  the  entire  three-dimensional  field. 

Recording  Holographic  Interferograms 

The  holographic  system  for  recording  inter¬ 
ferograms  near  a  model  rotor  blade  was  assembled 
in  the  Aeromechanics  l.aboratory  Anechoic  Hover 
Chamber.  Figure  5  shows  a  sehematte  of  the  opti¬ 
cal  3ystem  and  Fig.  6  shows  the  Anechoic  Hover 
Chamber .  A  ruby  laser  with  a  20-nsee  pulse  width, 
a  699.3-nm  wavelength,  and  a  power  of  1  J 
"freezes"  the  rotating  blade  at  any  desired  azi¬ 
muthal  angle.  A  beam-splitter  divides  the  laser 
beam  into  two  separate  beams  at  the  laser  outlet. 

A  microscope  objective  lens  expands  the  object 
beam  to  fill  a  2V- in . -dl am  spherical  mirror. 

.‘lirice  t tie  foci  of  both  the  objective  lens  and  the 
spherical  mirror  coincide,  a  collimated  plane  wave 
forms  as  the  beam  passes  through  the  rotor  area. 
The  object  beam  then  strikes  a  second  29-in. -diam 
spherical  mirror,  emerges  a3  a  converging  wave, 
and  illuminates  a  9-in.  by  b-in.  photographic 
plate.  The  reference  beam  is  lengthened  by  caus¬ 
ing  it  to  strike  several  plane  mirrors.  This  beam 
must  be  lengthened  so  that  the  dirrerence  in  the 
path  lengths  of  the  object  and  reference  beams  is 
1 033  than  the  coherence  length  of  the  laser  (one 
of  the  very  few,  and  easily  met,  geometrical  con¬ 
straints  in  a  holographic  system).  The  reference 
beam  is  expanded  by  an  objective  lens,  then  colli¬ 
mated  with  a  b-in.-diam  lens;  finally,  it  is 
directed  toward  the  film  so  that  it  overlaps  the 
object  beam. 

The  entire  procedure  is  conducted  from  out¬ 
side  the  hover  chamber,  once  the  optical  system  Is 
aligned.  Firing  the  laser,  changing  the  photo¬ 
graphic  film  plates,  and  controlling  the  test 
conditions  are  all  done  by  remote  control .  Recall 
that  to  record  an  interferogram,  two  exposures  at 
different  times  (different  flow  states)  must  be 
made  on  a  single  film  plate.  The  film  records  the 
first  exposure  while  the  rotor  blade  remains  sta¬ 
tionary.  In  this  case,  the  air  has  no  velocity 
and  therefore  ha3  a  uniform  refractive-index  dis¬ 
tribution.  The  film  records  the  second  exposure 
while  the  blade  rotates  at  the  desired  speed.  The 
nonhomogeneous  rerractive-index  distribution  in 
this  case  introduces  phase  changes  in  thp  second 
object  wave,  producing  Interference  fringes  on  the 
film  plate.  Thi3  double-exposure  recording  proce¬ 
dure  repeats  at  various  angles  around  the  flow  by 
synchronizing  the  laser  pulse  with  the  desired 
blade  position.  Because  of  the  long  optical  pelh- 
lengths  (90  ft),  the  recording  system  is  very 


sensitive  to  vibrations  or  the  optical  components. 
At  several  azimuthal  angles,  it  was  necessary  to 
record  multiple  interferograms  to  obtain  one  high- 
quality  i nterf erogram .  The  photographic  plates 
are  then  removed  from  the  recording  system  in  the 
hover  chamber,  developed,  and  reconstructed  in  a 
laboratory  for  further  processing. 

Holographic  interferograms  record  the  flow 
near  a  hovering  1 /7-scale  (geometric)  model  UH-1H 
rotor  with  untwisted  NACA.001?  airfoil  sections. 
The  blade  runs  at  a  tip  Mach  number  of  0.90  so 
that  the  flow  is  transonic  and  a  shock  wave  is 
present.'  This  model  normally  uses  two  blades; 
however,  in  view3  along  the  span,  the  optical  beam 
would  pass  through  the  refractive-index  field  of 
both  blades.  Because  the  refractive-index  fields 
of  the  two  blades  are  inseparable  at  these  angles, 
a  single-bladed  rotor  with  a  counterbalance  is 
used  instead  (Fig.  7). 

Holographic  interferograms  near  a  transonic 
rotor  blade  are  recorded  at  U0  different  viewing 
angles.  The  blade  rotates  in  a  clockwise  direc¬ 
tion  and  can  be  captured  at  any  desired  viewing 
angle  with  the  pulsed  laser.  The  tomography  code 
requires  flow  data  from  certain  viewing  angles 
within  a  180°  range.  Numerical  simulation 
results2^  using  numerical  computations  of  the 
riow3  suggest  recording  interferograms  from 
0  ■  8”  to  0  *  90°  and  from  0  *  190°  to 
n  =■  186°  in  2°  Increments,  as  defined  In  Fig.  8. 
The  missing  views,  8  *  92°  to  8  -  138°,  were 
presumed  to  have  very  few  interference  fringes  and 
were  not  utilized. 

Illustrated  in  Fig.  8  are  examples  of  holo¬ 
graphic  i nterf erograms  recorded  near  the  model 
blade  tip.  The  fringe  pattern's  appearance 
depends  on  the  viewing  direction.  Interferograms 
recorded  along  the  chord  (near  0  »  90°)  display 
very  few  interference  fringes,  since  the  optical 
rays  pass  through  the  field's  thinnest  (weakest) 
region.  No  observable  details  are  present  in 
these  views.  However,  in  views  along  the  span 
(near  0*0”  or  0  *  180“),  numerous  fringes  are 
visible,  because  the  optical  rays  pass  through  the 
longest  (strongest)  region  within  the  field.  The 
leading-edge  stagnation  point,  shock  structures, 
boundary-layer  separation,  ami  wake  system  are 
clearly  visible.  In  particular,  a  lambda  shock 
( il  *  180°)  and  the  radiated  shock  (0  *  186”) 
appear  above  the  blade.  Several  interferograms 
are  described  in  detail  in  Reference  2U. 

Pat  a  Extraction 

One  important  step  between  recording  the 
interferograms  and  applying  tomographic  recon¬ 
struction  techniques  is  the  evaluation  of 
the  interferograms;  that  is,  reading  fringe 


positions  and  fringe  numbers.  Up  to  now,  most 
I nterf erograms  from  aerodynamic  tests  have  been 
evaluated  by  either  reading  fringe  numbers  and 
their  positions  manually  or  by  tracing  the  fringe 
lines  by  hand  with  tire  help  of  a  tracking  device 
(for  instance  a  graphic  tablet).  Manual  evalua¬ 
tion  is  a  very  time-consuming  and  inaccurate  pro¬ 
cedure,  however,  and  it  is  evident  that  In  the 
current  application,  where  large  numbers  of  inter- 
ferograms  have  to  be  evaluated  at  several  hori¬ 
zontal  planes,  utilization  of  an  automatic  fringe¬ 
reading  procedure  would  enhance  the  evaluation  and 
would  make  the  interferometric  technique  a  much 
more  powerful  measurement  tool  (see  Reference  22 
for  a  detailed  discussion  of  this  procedure). 

In  order  to  reconstruct  the  three-dimensional 
flow  field  from  these  interferograms,  using  tomo¬ 
graphic  reconstruction  techniques,  the  fringe 
number-functions  along  cross  sections  of  a  plane 
parallel  to  the  rotor  disk  are  required.  Data 
f rom  all  the  views  are  needed  to  reconstruct  the 
index  of  refraction  in  a  particular  plane  above 
the  blade.  Several  planes  have  to  be  recon¬ 
structed  to  get  a  three-dimensional  flow-field 
representation.  This  procedure  requires  that  each 
i nterferogram  be  represented  in  a  form  that  allows 
the  fringe-order  function  to  be  computed  at  any 
desired  point. 

The  evaluation  of  i nterf erograms  using  com¬ 
pute; — aided  methods  can  be  subdivided  into  the 
following  steps:  1)  digitization  of  the  inter- 
ferograms  and  image  enhancement,  2)  fringe  seg¬ 
mentation  and  fringe-coordinate  extraction, 

3)  merging  of  fringe  fields  obtained  from  several 
magnified  views,  M)  fringe  numbering  and  correc¬ 
tion  of  fringe  disconnections,  5)  coordinate 
transformations,  6)  interpolation  and  extrapola¬ 
tion  of  fringe-number  functions. 

Hardware  Components.  An  image-processing 
system  (De  Anza  IP-6M00)  connected  to  a  VAX  11/780 
host  computer  provides  the  main  hardware  necessary 
for  digitizing  the  interferograms  and  for  doing 
some  image-enhancement  processing. 

The  resolution  or  the  system  is  51?  «  51? 
pixels  with  an  8-bit  intensity  range.  Currently 
it  Is  equipped  with  two  of  possible  four  memory 
planes,  as  well  as  with  a  graphic  and  an  alpha¬ 
numeric  overlay.  A  f rame -grabbi ng  unit  can  digi¬ 
tize  a  frame  of  a  video  signal  In  real  time 
(1/30  sec).  A  black  and  white  video  camera  (MTI 
series  88)  with  a  resolution  (bandwidth)  of  18  MHz 
is  connected  to  this  input  channel. 

The  system  ha3  an  arithmetic-logic  unit. 

(ALU),  with  which  real-time  addition,  subtraction, 
or  comparison  of  one  or  more  image  planes  may  tie 
made.  The  contents  of  each  memory  plane  may  be 


5 


routed  through  lookup  tables  before  being  input  to 
the  ALU,  to  the  video-output  processor,  or  to 
another  plane.  The  actual  contents  of  any  image 
plane  or  of  a  combination  of  image  plane3  is 
output  via  a  video-output  processor  and  can  be 
shown  on  a  color-display,  each  channel  has  its 
own  color-mapping  tables.  A  joystick  control 
device  is  used  for  interactive  input.  It  controls 
two  cursors,  which  may  be  used  in  a  number  of 
different  operating  modes.  A  color  print  system 
(Durm  Instruments  model  6J1  )  serves  as  a  hardcopy 
device  for  the  color  monitor. 

Digitization  and  Processing.  During  the 
recording  of  the  holographic  i nterf erograms ,  two 
fiducial  points  were  marked  in  the  image  plane 
willi  a  known  position  relative  to  the  rotor 
system.  This  allowed  a  subsequent  coordinate 
transformation  into  a  space-fixed  coordinate 
system.  Each  i nterf erogram  was  then  digitized, 
and  several  subsections  were  enlarged,  depending 
on  the  fringe  density  in  the  i nterf erograms .  A 
scale,  aligned  to  the  fiducial  marks,  was  always 
digitized  together  with  the  i nterferogram  so  that 
the  position  and  magnification  of  the  enlarged 
segments  could  be  identified.  Before  segmentation 
and  fringe-coordinate  extraction  are  applied,  it 
may  be  useful  to  do  some  image  enhancement  to 
reduce  noise.  Frame-averaging  methods  increase 
the  si gnal -to-noise  ratio  of  the  imaging  elec¬ 
tronic  components,  thereby  improving  the  picture 
qual i ty . 

One-Dlmensl onai _ Evaluation.  A  one- 

dimensional  interactive  processing  may  be 
employed.  If  only  one  two-dimensional  plane  of  the 
flow  field  is  reconstructed;  hence,  only  one  cross 
section  through  each  [nterferogram  has  to  be 
known.  In  this  case  some  of  the  expense  required 
for  the  full  two-dimensional  evaluation  of  the 
fringe  pattern  may  be  avoided.  Rut  the  procedure 
has  to  be  inter  ictive,  because  local  information 
(i.e.,  knowledge  of  the  fringe  locations  along  ore 
line)  is  not  sufficient  to  number  the  fringes 
correctly  or  to  detect  lost  fringes  or  other*  dis¬ 
tortions.  In  an  interactive  procedure,  informa¬ 
tion  such  as  acceptance  of  the  segmentation  and 
assignment  of  fringe  numbers  has  to  be  given  by 
the  user. 

A  program  featuring  t.li"  one-dimensional  eval¬ 
uation  digitizes  and  preproces.es  an  i  nterferogram 
and  does  a  fringe-segmentation  along  .a  line  or 
along  a  set  of  lines  through  the  field  of  view. 

The  result  of  the  segmentation  procedure  (a  binary 
fringe  pattern)  is  written  back  to  the  image 
screen  for  monitoring  reasons  (see  Figs.  9a 
and  9b).  The  user  may  interactively  change  the 
acceptance  threshold  or  edit  the  segmented  cross 
section  in  heavily  distorted  regions  before  he 
continues  to  the  numbering  section.  The  cursor 


can  be  moved  to  each  part  of  the  segmented  line, 
and  fringe  numbers  may  be  assigned  interactively 
by  various  commands.  The  task  of  assigning  fringe 
numbers  is  supported  by  color  coding  the  black 
parts  of  the  fringes;  this  shows  whether  the 
fringe-order  function  increases  or  decreases  by  1, 
or  if  a  discontinuity  is  present  between  adjacent 
fringes.  The  resulting  fringe-order  function  may 
also  be  plotted  onto  a  graphic  terminal.  A  post¬ 
processing  program  merges  all  the  data  taken  from 
different  enlarged  portions  of  an  i nterferogram 
(a3  shown  in  Figs.  9a  and  9b)  and  uses  a  spline- 
approximation  to  interpolate  between  the  fringes. 
The  final  result  is  extracting  the  fringe  order 
function  at  a  desired  height  above  the  blade  sur¬ 
face  as  shown  in  Fig.  9c. 

Two-Dimensional  Interactive  Processing. 

Often  the  knowledge  of  only  a  cross  section  across 
the  field  of  view  Is  not  sufficient  for  interpret¬ 
ing  a  given  flow  problem;  instead,  the  fringe 
patterns  have  to  he  evaluated  oven  the  whole  field 
of  view.  To  facilitate  the  two-dimensional  inter- 
ferogram  evaluation,  additional  algorithms  have 
been  implemented  to  segment  fringes,  to  extract 
each  fringe's  coordinates  (trace  the  rringe 
sides),  and  to  represent  the  whole  fringe  field  as 
a  polygonal  data  structure.  Methods  for  numbering 
these  line  fields  and  detecting  extraction  errors 
have  been  developed.  With  a  two-dimensional 
Interpolation  scheme,  fractional  fringe-order 
numbers  may  then  be  estimated  at  any  given  point 
in  the  field  of  view. 

Figure  10  shows  an  example  of  a  typical 
fringe-coordinate  extraction.  An  i nterferogram 
(Fig.  10a)  Is  digitized  over  the  whole  field  of 
view  with  low  resolution.  The  white  lines  over¬ 
laying  the  display  represent  the  extracted  fringe 
sides,  and  the  black  line  represents  the  border 
polygon.  The  border  polygon  was  input  by  hand 
before  the  coordinate  extraction  in  order  to  mask 
off  background  objects  and  the  part  of  the  fringe 
field  with  high  fringe  density.  (In  our  implemen¬ 
tation,  all  the  graphics  are  in  color  and  there¬ 
fore  easier  to  distinguish  from  the  displayed 
image. ) 

As  was  discussed  before,  it  may  be  necessary 
to  digitize  enlarged  subsections  to  resolve 
regions  with  high  fringe  densities.  The  resulting 
fringe  polygonal  field  is  then  combined  out  of  all 
the  individual  digitized  and  processed  subsec¬ 
tions,  starting  with  the  one  with  highest  resol u 
tion.  To  handle  the  problem  of  partly  overlapping 
polygonal  line  fields,  a  boundary  polygon  is  main¬ 
tained  tnr  each  subsection.  This  defines  the 
definition  area  for  each  fringe  polygonal  field. 
Upon  combination  of  two  overlapping  rields,  the 
border  polygon  or  the  one  with  higher  resolution 
(priority)  is  used  to  intersect  the  lines  of  the 


other  field.  Both  line  fields  are  then  connected 
at  those  Intersection  points.  The  new  border 
polygon  i3  the  border  of  the  combined  areas. 

Figure  10b  shows  an  enlarged  digitized  part  of 
the  Interferograra  shown  in  Fig.  10a,  with  the 
extracted  fringe  polygons  overlaying. 

CAT  Reconstruction 

Fringe-order  functions  are  transferred  to  the 
filtered  back-projection  CAT  code  which  computes 
the  refract  1 ve-lndex  field  at  specific  points  in  a 
chosen  horizontal  plane  above  the  blade  surface. 
The  code  assumes  refractionless  light  rays;  there¬ 
fore,  each  horizontal  plane  can  be  treated  inde¬ 
pendently,  even  though  data  for  each  plane  is 
taken  from  one  set  of  interferograms.  The  pertur¬ 
bation  velocity  is  computed  by  first  converting 
refractive  index  to  density,  using  the  Gladston- 
Dale  relation: 

p  »  (n  -  1 )/k 

Density  13  then  converted  to  perturbation  velocity 
from  a  form  of  Bernoulli's  equation  for  steady 
(with  respect  to  the  rotation  blade),  compressi¬ 
ble,  l3entropic  flow: 

V  -  -OR  *  |(«R)2  -  y~~T  [(^)  *  ]|>  <*> 

The  procedure  is  repeated  in  several  planes  above 
the  blade  to  reconstruct  the  entire  three- 
dimensional  field  near  the  model  blade  tip. 


Reconstruction  Results 

The  holography-CAT  reconstruction  of  the 
blade-t' p  velocity  field  is  compared  with  numeri¬ 
cal  computations.  The  computations  used  here  are 
conservative,  mixed-difference  solutions  of  the 
transonic  smal 1 -disturbance  equation.  Results 
from  both  sources  are  presented  in  four  horizontal 
planes,  as  Identified  In  Fig.  11a.  Three  plot 
types  are  used  to  visualize  the  flow  field. 

First,  velocity  contours  are  given  In  plan  view 
(see  Fig.  11b),  where  the  blade's  leading  and 
trailing  edges  are  at  X/C  *  0.0  and  X/C  *  1.0, 
respectively.  The  blade  tip  is  located  at 
R/R(j  *  1.0,  the  rotation  center  is  at  R/Rq  ’  0-0, 
and  the  blade  rotates  in  a  clockwise  direction. 
Oeeond,  perspective  views  are  displayed  in  which 
velocity  val"°3  are  plotted  along  the  vertical 
axis.  The  d  ita  and  geometry  are  identical  to  the 
contour  plots,  though  the  data  are  viewed  from 
near  the  rotor  hub.  Third,  velocity  distributions 
are  shown  at  six  radial  locations  (see  Fig.  lib) 
for  each  plane. 


Figure  12  compares  the  velocity  contours 
derived  from  the  holograph! c-CAT  method  and  the 
numerical  computations  near  the  blade  surface 
(Y/C  -  0.08).  Both  methods  display  low-velocity 
regions  near  the  leading  and  trailing  edges,  and 
also  display  a  high-velocity  region  over  the  blade 
surface  containing  a  shock  at  approximately 
X/C  -  0.60  near  the  blade  tip.  The  general  con¬ 
tour  shapes  show  a  strong  resemblance  except  near 
the  blade  tip  (roughly  the  last  61  of  blade  span). 
The  maximum  velocity  region  appears  at  the  shock 
foot  In  the  computational  analysis,  but  It  appears 
closer  to  the  leading  edge  and  farther  from  the 
shock  foot  in  the  holography-CAT  results. 

Figure  13  shows  velocity  values  for  the  same 
plane  in  perspective  view.  Again,  the  general 
flow  shapes  appear  very  similar.  The  major  dif¬ 
ference  between  the  two  results  is  the  roughness 
(minor  ’’ridges")  in  the  reconstructed  flow.  This 
may  be  due  to  reconstruction  artifacts  caused  by 
nonconti nuous  data  (interferograms  recorded  in  ?° 
azimuthal  increments)  or  by  noise  (erroneous 
fringes)  in  the  lnterferogram  data  caused  by  opti¬ 
cal  component  motion.  Figure  14  compares  velocity 
distributions  at  six  radial  locations.  The  rough¬ 
ness  of  the  reconstructed  now  can  be  seen 
throughout  the  figures.  Also,  the  holography-CAT 
method  determines  the  shock  location  to  be 
slightly  more  downstream  (3t)  than  does  the  numer¬ 
ical  code  solution.  The  major  differences  ran  be 
observed  in  Fig.  I4d,  where  the  discrepancies  at 
the  leading  edge  and  over  the  blade  surface  are 
clearly  visible.  The  leading  edge  (X/C  *  0.0) 
difference  may  result  from  a  breakdown  of  the 
small  crossflow  assumption  [Eq.  (4)]  near  the 
blade  tip.  The  difference  over  the  blade  surface 
(X/C  -  0.2  to  X/C  *  0.6)  may  be  attributed  to 
the  existence  of  a  lambda  shock,  in  the  lnterfero- 
gram  data  (t.e.,  Fig.  8f ) ,  which  cannot  be  pre¬ 
dicted  by  the  nonviscous  numerical  potential  code. 
A  lambda  shock  was  also  observed  in  Schlleren 
photographs  from  a  previous  wind-tunnel  test 
(Fig.  16,  taken  from  Reference  26)  using  the  same 
airfoil  and  tip  Mach  number. 

Figure  16  compares  the  velocity  contours  Trom 
both  the  holography-CAT  method  and  numerical  code 
at  Y/C  -  0.22  (near  the  upper  region  of  the 
lambda  shock).  The  general  velocity  contours  show 
an  excellent  agreement  In  both  shape  and  magnitude 
throughout  the  plane.  The  maximum  velocity  region 
on  the  blade  surface  (near  X/C  *  0.60, 

R/Rq  -  0.96)  match  much  closer  in  this  rlane  than 
in  the  plane  near  the  blade  surface  (Fig.  12). 

The  perspective  view  of  Fig.  17  also  shows  an 
excellent  agreement  in  shape  and  magnitude. 

Again,  the  most  noticeable  difference  is  the  extra 
ridges  in  the  reconstructed  flow  at  the  name  loca¬ 
tions  and  orientations  as  seen  in  the  previous 
plane  (Fig.  13a).  The  velocity  distributions  in 
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this  plane  (Fig-  18)  compare  favorably,  especially 
inboard  of  the  blade  tip  (see  Figs.  I8a-l8cl,  as 
well  as  at  the  blade  tip  (Fig.  l8d).  The  lambda 
shock's  effect  near  the  tip  is  apparently  weaker 
in  this  plane;  thus,  the  velocity  distribution 
magnitudes  and  shapes  are  much  closer  than  those 
i  n  F  i  g .  1  *4 . 

Figure  Id  compares  velocity  contour's  at 
Y/C  =  0.49  above  the  blade.  The  velocity  contour 
shapes  and  magnitudes  are  similar  in  all  regions 
except  that  the  hoi ography-CAT  method  shows  the 
maximum  velocity  point  to  be  slightly  (l£  span) 
outboard  of  the  numerical  result.  The  extra 
ridges  soon  in  the  perspective  view  (Fig.  20)  are 
in  approximately  the  same  location  and  orientation 
as  those  in  the  previous  plane;  however,  the  mag¬ 
nitude  of  the  ridges  has  decreased.  Figure  21 
•.Mows  velocity  distributions  for  this  plane. 

There  is  good  agreement  at  all  locations  except 
beyond  the  blade  tip  (Fig.  21 e)  where  the  hoi og¬ 
raphy-CAT  method  shows  slightly  larger  velocity 
rnagn i tudes . 

Finally,  Fig.  ??  shows  velocity  contours 
for  Y/C  1.17.  The  velocity  contours  match 
throughout.,  though  there  are  no  distinguishing 
features  in  the  flow  at  this  height  above  the 
hi  ride  surface.  In  this  plane,  the  ext  ra  ridges 
arc  almost  unobservable  in  the  perspective  view 
(Fig.  21),  The  velocity  distributions  of  Fig.  24 
also  show  a  strong  .similarity  between  the  two 
methods  at  all  radial  locations,  for  both  results 
show  that,  this  plane  is  at  the  perturbed  flow's 
ip  per  extent. 

Overall,  the  agreement  between  the  holog- 
raphy-t’AT  results  and  the  numerical  solution  is 
•mV nemo 1 y  eneouragi ng .  However,  comparisons  with 
other  experimental  data  sources  are  required 
before  a  final  decision  can  be  made  about  the 
ho] ography-CAT  results.  Several  discrepancies 
must  be  res  Ived.  First,  data  from  prossure- 
i  nst.rumcnted  blades  will  aid  in  confirming  the 
shock  location  and  whether  a  lambda  shock  does 
#'xist.  in  the  flow.  Second,  the  extra  ridges 
appearing  in  the  reconstruct ed  results  must  be 
eliminated  or  reduced;  the  ridges  may  be  caused  by 
poor  data  recorded  in  the  l nterf erograms .  To 
Improve  the  quality  of  the  I nterferograms ,  a  modi¬ 
fied  optical  system  (for  both  hover  and  forward 
flight  testing)  is  necessary.  Tn  addition, 
neither  the  h>>l ography-CAT  reconst ruct l on  nor  the 
fiumcrie.il  code  solut  ions  compute  the  expected 
shock  strength  beyond  the  blade  tip.  Acoustic 
measurements  indicate  a  much  stronger  radiated 
shock  than  is  indicated  by  these  two  results. 


Concluding  Remarks 

The  holographic  i nterferometry  computer- 
assisted  tomography  technique  proved  to  be  a 
highly  effective  way  of  measuring  the  three- 
dimensional,  transonic  flow  field  near  a  model 
rotor-blade  tip.  Results  from  this  method  compare 
favorably  with  those  of  numerical  computations, 
except  very  near  the  1. 1  p  region.  That  discrepancy 
may  be  due  to  the  existence  of  a  lambda  shock 
recorded  by  the  i nterferograms  (and  In  previous 
Schlleren  photographs)  which  is  not  predicted  by 
the  nonvi  scour,  potential  code.  In  other  regions, 
the  velocity  distributions  along  the  chord  are 
similar  in  both  shape  and  magnitude.  However,  the 
results  from  the  technique  must  he  further  verl 
fled  against  other  experimental  data. 

Since  this  is  the  first  successful  Implemen¬ 
tation  of  the  holograph! e  i nterf eromet ry ,  com¬ 
puter-assisted  tomography  method  in  rotor  flow 
studies,  many  improvements  are  indicated.  For 
example,  t.he  optical  system  must  be  improved  so 
that  better  quality  i  nterferograms  can  lie 
recorded,  and  an  automatic  f r 1 nge-readi ng  tech¬ 
nique  must  be  completed  so  that  the  time  required 
to  evaluate  i nterferograms  can  be  shortened.  Upon 
verification  of  these  results  and  after  the  system 
is  improved,  measurements  of  other  model  rotor- 
blade  flow  fields,  including  those  of  forward 
flight,  will  bo  performed. 


Append! x  : _ Holograph i c  I nterf prompt ry 

Doubl e- exposure  holographic  i nterf erometry  Is 
the  i nterf eromotr i c  comparison  of  two  object  waves 
that  are  recorded  holographically  sequentially  in 
time.  The  i nterf erogram  is  recorded  by  first 
exposing  a  photographic  film  to  a  reference  wave 
and  an  "undisturbed"  object  wave,  as  shown  In 
Fig.  25a.  batep  in  time,  the  same  photographic 
plate  is  exposed  to  a  reference  beam  and  a  second 
"disturbed"  object  wav.,  as  shown  in  Fig.  ?Sb . 

The  irradiancc  at  the  film  plane  is 
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The  last  two  terms  of  Fq .  ( A 1 )  represent  the 
i  nterf  er’ence  pattern  recorded  on  the  film,  which 
contains  both  amplitude  and  phase  information 
about  the  two  reference  and  t.he  two  object  waves. 
The  amplitude  transmittance  of  the  developed  film, 
called  a  holograph. c  l nterf erogram ,  is 
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When  the  reconstruction  wave  illuminates  the 
i nterferogram,  as  shown  In  Fig.  25c,  the  trans¬ 
mitted  light  is 

Ui  '  Ucl 

'  V2Ir  +  !01  *  V  *  WU01  +  V 
*  UrVU01  *  V‘  (A2) 

The  second  term  in  Eq.  {A 2)  is  the  "composite" 
virtual  image  reconstruction  of  the  two  object 
waves  (Uqj  *  Uq^ ) •  By  illuminating  the  inter- 
ferogram  with  a  conjugate  reconstruction  wave 
(Fig.  25d),  the  transmitted  light  is 


numerical  data  shown  in  this  paper.  We  thank 
Professor  L.  He3selink  of  Stanford  University  for 
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guidance. 
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Fig.  15  Schlieren  photographs  illustrating  lambda  shock  on  a  NACA  0012  airfoil  operating  at  a  tip  Mach 

number  of  0.90  in  forward  flight.  (From  Ref.  25.) 


Fig.  16  Perturbation  velocity  contours  In  plan  view  for  Y/C  =  0.22  above  blade  centerline: 

interval  ■  50  ft/sec.  a)  Holography-CAT  reconstruction;  b)  numerical  code. 


contour 
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Fig.  12  Perturbation  velocity  contours  in  plan  view  Tor  Y/C  -  0.08  above  blade  centerline:  contour 
interval  -  50  ft/sec.  a)  Holography-CAT  reconstruction;  b)  numerical  code. 


<«)  (b) 

Fig.  13  Perturbation  velocity  values  Tor  Y/C  -  0.08  above  blade  centerline,  a)  Holography-CAT 

reconstruction;  b)  numerical  solution. 
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Fig.  6  Holographic  setup  at  Anechoic  Hover 
Chamber. 


Fig.  7  One-blade  rotor  with  a  counterweight 
balance. 


Fig.  3  Filtered  back  projection,  a)  The  projection  data  are  convolved  (filtered)  with  a  suitable  process¬ 
ing  function  before  back-projection;  b)  three  back-projections  of  an  absorbing  disk.  (From  Ref.  19). 


Fig.  4  Recording  interferograras  at  various  angles  around  the  field  of  interest  for  tomographic 

reconstruction. 


Fig.  25  Double-exposure  holographic  interferometry .  a)  First-exposure  recording;  b)  3econd-expos 
recording;  c)  reconstruction  of  the  virtual  Image;  d)  reconstruction  of  the  real  image. 
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